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R: 


The  approximate  theory  of  radar  reflection  from  thin  wires  by 
Van  Vleck,  et.  til. ,  gives  very  good  average  radar  cross  section  (RCS) 
results  and  good  angular  RCS  results  except  end -on.  In  this  paper,  the 
nature  of  this  end -on  discrepancy  is  examined.  It  is  found  that,  if  the 
complete  expressions  derived  by  Van  Vleck,  et.  al. ,  are  ut&ixed  without 
the  approximations  made  to  simplify  calculation  of  the  average  RC5<  over  all 
angles  of  incidence,  then  very  accurate  RCS  results  are  predicted  for  nearly 
all  angles  of  incidence  and  for  all  the  wire  length -to -wavelength  ratios  between 
0.45  and  50.  Comparisons  with  the  numerical  results  of  a  source  distribution 
technique  (SDT)  computer  program  and  with  the  results  due  to  Ufimtsev  are 
shown. 
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I.  INTRODUCTION 


In  a  previous  report  (Ref.  1),  predictions  of  the  radar  cross 
section  (RCS)  of  long,  thin  wires  by  several  authors  were  compared  with 
experimental  data.  It  was  found  that  the  predictions  of  Van  Vleck,  et  al. , 
(Ref.  2)  were  of  greater  validity  than. had  previously  been  believed;  however, 
the  RCS  of  thin  wires  for  end-on  incidence  was  considerably  in  error.  This 
error  was  believed  to  be  due  to  the  approximations  made  by  Van  Vleck,  et  al. , 
to  simplify  the  calculation  of  the  average  RCS  over  all  angles  of  incidence. 

To  confirm  this  belief,  the  theory  of  Van  Vleck,  et  al. ,  is  examined 
in  detail  in  Sec.  2  of  this  report.  It  is  shewn  that,  without  some  of  these 
simplifying  approximations  discussed  above,  the  general  theory  derived  by 
Van  Vleck,  et  al. ,  does  indeed  give  a  cross  section  that  goes  to  zero  at  end- 
on  incidence.  The  general  theory  is  used  to  calculate  the  RCS  of  thin  wires 
of  various  lengths.  These  results  are  compared  with  data  generated  by 
BRACT  ,  a  source  distribution  technique  (SDT)  computer  program,  to  verify 
the  fact  that  use  of  the  general  theory  instead  of  the  approximate  theory 
(Ref.  2)  has  not  degraded  the  overall  angular  RCS  results. 

The  approximate  theory  of  Van  Vleck,  et  al.  ,  is  also  compared 
directly  with  the  general  theory  and  with  BRACT  to  determine  the  differences 
among  the  results.  Note  that  resalts  obtained  using  the  theory  of  Ufimtsev 
(Refs.  1  and  3)  have  also  been  compared  with  the  BRACT  results.  These 
comparisons  are  presented  in  Sec.  3  of  this  report.  A  brief  discussion 
follows  (Sec.  4). 


The  BRACT  computer  program,  which  solves  'V  >  thin  wire  integral  equation 
of  the  complete  electromagnetic  scattering  problem  for  arbitrary  wire 
structures,  hai.  been  validated  through  extensive  >’se  and  is  considered  <.;> 
have  an  accuracy  of  better  than  1  IB  (Refs.  1  and  4). 


II.  THEORY 


Consider  a  plane  wave  incident  on  a  thin  wire  of  length  21,  as 
shown  in  Fig.  i . 


Fig.  1.  Plane  of  Incidence 


The  angle  of  incidence  is  6,  and  the  angle  between  Eq  (the  electric  field 
vector}  and  the  plane  formed  by  the  wire  and  the  propagation  vector  p  is  0. 
The  CGS  system  of  units  is  used  throughout  this  paper.  Then,  for  the 
scattered  field  (EL),  from  Maxwell’s  homogeneous  equations 


i  3A9 
c  at 


(t) 


■ . »  g  g 

where  A  is  the  scattered  vector  potential  and  <p  is  the  scattered  scalar 
potential.  If  the  Lorentz  gauge  is  used 


V"*  A  +  - 

c 


(2) 


Ffgesding  mm 
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Eq.  (1)  becomes 


ii^  =  cvv.xs.i  s-% 

dt  CV  V  A  c  ^ 


(3) 


If  all  time  dependence  appears  at  a  single  frequency  as  e”lwt,  then  Eq.  (3) 
becomes 


■iwE  =  cvr  •  A  +  ~A 

c 


(4) 


Now,  if  only  the  tangential  component  of  the  field  at  the  surface  of  the  wire 

is  considered  and  if  it  is  assumed  that  the  wire  is  a  perfect  conductor  lying 

on  the  z  axis,  E8  =  -  E1  (E  1  is  the  incident  field)  and  Eq.  (4)  becomes 
z  z 


82A  - 

tt*  +  e  \ 

a  z 


(5) 


where  the  superscript  s  has  been  dropped.  The  incident  field  is  given  by 


E  e 
o 


ip  •  X  -  iwt 


If  the  time  component  is  neglected  (Fig.  i) 


E1  =  E  cos  O  sin  8  e 
z  o 


ipz  cos  0 


(6) 


-4- 


or 


92A 


+  p^A  =  ipE  cos  0  sin  0  e***z  cos  ® 

a  z  2  ° 


The  homogeneous  solution  of  Eq.  is 


A  cos  pz  +  B  sin  pz 


and  the  inhomogeneous  solution  is 


iEQ  cos  0  sin  oj"  e*^  COS  ®  sin  p(z  -  |)  d£ 


i£  cos  0 
o 

p  sin  0 


(e1^2  -  cos  pz  -  i  cos  G  sin  Pz) 


where  q  =  p  cos  G  =  2it/\  cos  G. 

Recall  from  Maxwell's  inhomogeneous  equations  that  the 
potential  due  to  a  cnrrent  distribution  is  given  by 


v2a  *  ^ 

c  dt 


4tt 


which  has  a  general  free  space  solution 


A  =  i  -  x;|/c)  t]  dti  dx,3 

!  X  ~  X  '  I 


(7) 


(8) 


(9) 


vector 


(iO) 


(H) 


If  J  =  J(X)e”lwt  and  if  it  is  assumed  that  the  current  flows  only  along  the  z' 
direction  and  resides  at  the  axis  of  the  wire,  while  the  vector  potential  is 
wanted  at  the  surface  of  the  wire 

2  r  -iPr  / - 5 - J 

A  (a)  =  j  —————  dz'  and  r  =  *J(z  -  z')  +  a  (12) 

*  c  r 

2 

Or,  if  J  =  I /it a  and  if  the  vector  potential  solutions  of  Eqs.  (7)  and  (12)  are 
equated 


i&E  cos  $ 

dz‘  =  A.  cos  Ps  +  B.  sin  pz  + - =p - 

P*  sine 


eiqz 


(13) 


This  is  identical  to  Eq.  (1)  in  Ref.  2;  note  that  the  coefficients  of  the  sine 
and  cosine  terms  of  Eqs.  (8)  and  (9)  have  been  combined.  Equation  (13)  is 
an  integral  equation  that  must  be  solved  for  the  current  on  the  wire,  subject 
to  the  boundary  condition  that  the  current  vanish  at  ski.  Further,  the 
current  should  vanish  for  6  =  0,  i.e. ,  the  right-hand  side  of  Eq.  (13)  should 
vanish  at  6  =  0. 

In  order  to  solve  for  Eq.  (13),  assume  that  the  current  I(z)  is 

given  by 


I(z) 


a  el<*z  +  Yj  cos  pz  +  iY^  sin  pa 


(14) 


where  a,  Y^,  and  Y-,  do  not  depend  on  z.  Note  that  Eq.  is  just  an 
approximation  to  the  current  and  that  a  more  accurate  expression  would 
isquire  an  iterative  solution  (to  be  discussed  belcvj. 


-6- 


The  value  of  a  is  determined  by  substituting  the  el<*^  term  of 
Eq.  (14)  into  Eq.  (13)  and  by  equating  with  the  el<*z  terms 


J.t 


-iSr  iwE  cos  <b 

dz'  =  — £ - 

sin  6 


The  integral  on  the  left-1  >nd  side  is  broken  into  three  terms  for  e  *sy 
evaluation  (Ref.  5) 


fgl!!  dz'  S  e*2  £2Lfil  dm,  + 

r  «/-!  r  JLl  r 


cos  (3r  dz' 

(16) 


iqz'  •  « 

e  ^  sm  Br 


L“> 


dz'  =  Z(z)  = 


-  cos  Sr) 


dz  (17) 


The  second  term  on  the  right-hand  side  of  Eq.  (17)  will  not  go  to  zero  for 
r  =  0;  therefore,  we  may  replace  r  by  |  z  -  z'|  for  this  term  and  for  the  two 
right-hand  terms  in  Eq.  (16).  By  a  change  of  variables  in  Eqs.  (16)  and 
(17)  and  by  use  of  the  sine  and  cosine  integrals 


/x  , 

1  ••  cos  t 

— t — 


/X  .  . 

Sift  at 
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Equations  (16)  and  (17)  are  evaluated  as 


Z(z)  =  log 


[(I  +  z)2  +  a2ji/2  +  (I  +  z) 
[(.«  -  z)2  +  a2]1/2-  (|  -  Z) 


-  Cin  |3(l  +  z)  -  Cin  (5(1  -  z) 


(18) 


K(aJ  ft1  e*ipr  d*‘ 

=  2  <s*qS{2Z(a)  +  2 Cin  p(l  -  a)  +  2Cin  p<*  +  a)  -  Cin(P  +  q)(l  -  a)  -  Cin(p  -  q){<  +  z) 
••  C>n(P  -  q)(i  -  a)  -  Cin(p  +  q ){i  +  a) 

-  i  [Si(P  +  q)(l  -  a)  +  Si(P  -  q)(i  +  a)  +  Si(P  -  q)(i  -  a)  +  Si(P  +  q)(i  +  a)]  } 

+  j  eiq*{Cin(P  +  q)«  -  a)  +  Cin(P  -  q)(l  +  a)  -  Cin(p  +  q )(*  +  a)  -  Cin(P  -  q)(*  -  a) 

+  ifSi(p  +  q)(l  -  a)+  Si(P  -  q)(l  +  z)  -  Si(p  +  q)(l  +  a)  -  Si(P  -  q)«  »  a)  j  } 


Since 


iw  cos  <p  E  e 
oK(z)  el<*Z  =  - 


lqz 


p  sinG 


(19) 


and  since  a  should  be  independent  of  a,  an  average  over  the  wire  is  taken 


Hz)  =  2  jlog(2i/a)  +  5og  2  +  a/2 1  -  Cin  2PI  -  ] 

!?(*)  =  2  log{2i/a  )  +  2  log  2  +  a/I  -  Cin  2(P  +  q)l  -  Cin  2{p  -  q )/  - 


-  i  | Si  2<P  +  q)l  +  Si  2{P  -  q)J  +  + 


(20) 
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or 


a  - 


iw  cos  0  E 
_ T  c 

Kp2  sine 


(21) 


Similarly,  Yj  and  Y^  are  related  to  Aj  and  by 


Aj  =  LY^  and  =  iLYg 


(22) 


where 


L  =  2  log(2f/s)  +  2  log  2  +  a/I  -  Cin  4(31  -  { 

.l(Si4  pi+g>iffi-«) 


(23) 


Note  that  when  q  goes  to  {3  or  when  0  goes  to  zero,  the  value  of  K 
approaches  the  value  of  L.  The  use  of  the  complete  expressions  for  K  and  L 
differs  from  Van  Vleck's  usage  in  that  he  uses  the  asymptotic  values  of  the 
Cin  and  Si  functions  for  large  argument.  The  usage  is  incorrect  when 
P  ±  q  -*  0. 

In  order  to  determine  and  Bj  or  Yj  and  Y^,  the  boundary  con¬ 
ditions  on  the  current  are  applied.  However,  if  they'  are  applied  directly  to 
Eq.  (14),  the  results  would  give  an  infinite  current  for  (3  =  ±mr/2;  n  =  1,2,3 
•  *  •  and  is  identically  zero  for  q  =  ±p,  Instead,  the  integral  in  Eq.  (13)  is 
broken  into  three  terms,  as  in  Eq.  (16) 

I(z)  Z(z)  -  r  *  [i(z')  -  I(z)]  cob  pr  dz'  +  i  f  r  *  I(z')  sin  pr  dz': 

J-l  J~l 


+  /  \  Eq  elqs  +  Aj  coc  pz  +  Bj  sin  pz 

\p£  sin  0  /  ° 


(24) 
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It  has  previously  been  mentioned  that  a  more  accurate  representation 
of  the  current  would  require  an  iterative  solution.  Equation  (24)  would  be  the 
basis  for  such  a  solution,  with  the  previous  solution  (e.g.,  Eq.  14)  for  the 
current  used  in  the  right-hand  side  of  Eq,  (24).  The  averaged  value  of  Z(z) 
would  be  used  in  performing  any  integrations  past  Ihe  first  iteration  to  the 
current.  Instead  of  iterating,  Eq.  (24)  is  forced  to  obey  the  boundary  con¬ 
ditions  1(H)  =  0  by  using  the  zero4*1  approximation  to  the  current  (Eq.  14)  in 
the  right-hand  side  ol  Eq.  (24)  to  determine  and  V^.  This  leads  to 
Eqs.  (25)  and  (26) 

Q  eiql  +  At  cos  pi  +  Bt  sin  pi  =  aD  eiqi  +  (Yj/2)(  E  elPi  +  F  e~lP<) 

+  (Y2/2)(e  eipi  -  F  e”lpi)  (25) 

Q  e'iqi  +  At  cos  pi  -  Bj  sin  pi  =  aG  eiqi  +  (y1/2)(e  elpi  +  F  e'lPi) 

-  (y2/2)(e  eipi  -  F  e"lPi)  (26) 

where 

i  w  cos  0  E 

Q  =  — - 2.  (27) 

p  sin  6 

and 

D  =  Cin  2P1  -  Cin  2(P  +  q)i  -  i  Si  2(p  +  q)i 
G  r  Cin  2pi  Cin  2(P  -  q)i  -  i  Si  2(P  -  q)i 
E  =  Cin  2pi  -  Cin  4 pi  -  i  Si  4pf 

F  =  Cin  2pi  (28) 


-10- 


From  Eq.  (22),  the  following  solutions  are  obtained  for  Vj  and  Y? 


Vi  = 


St, 


V*""P 


and  the  current  is 


1(a)  =  S  |e*qZ  -  Tj  cos  (3z  -  iT^  sin  pzj 


(29) 


where 


T  =  [2K  cos  qf  -  (D  eiq*  +  G  e"lqi)] 
1  [2L  cos  pi  -  (E  ei(^  +  F  e~ipf)j 


(30) 


X  -  [2  iK  sin  ql  -  (D  e*q*  -  G  e  *q-*)] 
2  [2  iL  sin  Pi  -  (E  eiPi  -  F  e"1^)] 


(31) 


The  current  thus  correctly  goes  to  zero  for  q  =  ±p. 

The  vector  potential  in  the  far  field  is  then  given  by 


A(X) s  —  I(z')eiqVdz'j 


(32) 


where  k  is  the  unit  vector  in  the  z  direction  and  q'  =  (3  cos  0',  with  0'  the 
received  angle. 


The  far  field  scattered  E  (from  Maxwell's 


-*  •  _*  *a  ij3r  r‘ 

E  =  |7X7XAS*  — -  sin  8  'j 


equation)  is  written  as 
I(z')  eiq'z'dz'  (33) 


Assume  that  the  monostatic  cross  section  is  desired  and  that  the  polarization 
angle  of  the  detector  is  the  same  as  that  of  the  transmitter.  The  scattered 
field  detected  is  then 


•E  e 


>0* 


E  = 


i_  co„2Ji^ii .  [t, + t  j 

prj^  j^cos  8  1  2  1  +  cos  8 

.  [T  ..  T  I  (P  USLlil 
lll  2  i  -  cos  8  J 


(34) 


The  RCS  of  the  wire  is  defined  as 


or  (8,  4>) 


4irrZE2 


\ 


2  4 


cos 
ttK  K* 


where  Ej  is  the  quantity  in  brackets  in  Eq.  (34). 


(35) 
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HI.  COMPARISON  OF  RESULTS 


The  expressions  given  in  the  previous  section  have  beer  programmed, 
and  calculations  have  been  performed  for  the  cases  listed  in  Table  t . 


Tabie  i .  Thin  Wire  Parameters 


— 

r 

Polarization* 

i_ _ 

Case  No. 

Wavelength,  m 

Transmitted 

Received 

pa 

Pi 

1 

1.00 

Linear 

Linear 

3.14  X  10'2 

1.415 

2 

0.227 

Circular 

Circular 

4.2  X  10‘3 
-2 

4.  44 

3 

1.00 

Circular 

Circular 

3.95  X  10 

17 

4 

0.69 

Linear 

Linear 

9.1  X  10-3 

34.6 

5 

0.02 

Linear 

Linear 

4.78  X  10"2 

157 

Circular  transmitted  and  received  RCS  are  6  dB  lower  than  linear  trans¬ 
mitted  and  received  RCS. 


When  these  calculations  are  compared  with  data  generated  by  the 
BRACT  computer  program  (Ref.  1)  for  the  same  cases  (Figs.  2  through  6), 
it  is  seen  that  the  general  Van  Vleck  calculations  are  nearly  identical  to 
those  of  BRACT.  There  is,  however,  a  difference  of  at  most  2  dB  in  the 
maximum  RCS  in  the  lobes  between  end-on  (0  deg)  and  broadside  (90  deg) 
for  the  pi  =  17  and  157  cases.  Further,  the  nulls  for  the  general  Van  Vleck 
theory  are  much  deeper  than  for  BRACT  for  the  larger  kl  values.  A  com¬ 
parison  of  the  general  and  approximate  Van  Vleck  theories  is  presented  in 
Figs.  7  through  11.  Here  it  is  seen  that,  except  for  end-on,  the  results  are 
almost  identical  to  those  of  Figs.  2  through  6.  Note  that  the  approximate 
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Fig..  2.  Generalised  Van  Vleck  versus  SDT  ftesuits:;  Monostatic  Cross 
Section  of  a  Dipole,  Linear  Polarisation 
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.  Generalized  Van  Vleck  versus  SDT  Results:'  Monostatic  Cross 
Section  of  a  Wire,  Circular  Polarization,  2t/\  -  1.4097 
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Fig.  4.  Generalized  Van  Vleck  versus  SDT  Results:,  Monostatic 
Section  of  a  Wire,  Circular  Polarization,  21 /\  5,4113 


Fig.  5.  Generalized  Van  Vleck  versus  SDT  Results:  Monostatic  Cross 
Section  of  a  Wire,  Linear  Polarization,  21 /\  -  11.0145 
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Original  Van  Vleck  versus  Generalised  Van  Vlock  Results: 
Cross  Section  of  a  Dipole,  Linear  Polarisation 
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Original  Van  Vle-ek  varsus  Generalis'd  Van  VU;ek  Results:  Monost 
Cross  Section  of  a  W  .'re,  ‘  ircnlar  l-olnri*ntn:n,  Z  * /"'<  S.  41!  > 


Van  Vleck  theory  agrees  better  with  BRACT  (Fig.  12)  than  does  the  general 
theory  except  for  end- on  incidence  for  larger  pi  values.  This  in  seen  by 
comparing  Figs.  6  and  12.  For  the  purpose  of  comparison,  Ufimtsev's 
equations  have  been  programmed  and  used  to  calculate  the  RCS  of  die  thin 
wire  cases  given  in  Table  1 .  These  calculations  are  presented,  with  die 
BRACT  calculations,  in  Figs.  13  through  17.  It  is  seen  that  the  agreement 
with  BRACT  is  not  as  good  as  the  agreement  between  the  general  Van  Vleck 
and  the  BRACT  results  for  pi  <  35;  for  these  larger  pi  values,  results 
obtained  using  Ufimtsev’s  equations  generally  agree  more  closely  with 
BRACT  than  do  the  general  theory  results,  especially  in  die  RCS  nulls. 
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Fig.  13.  TJfitntsev  versus  SDT  Results:,  Monostalic  Cross  Section  of  a 
Dipole,  Linear  Polarization 
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Fig,  14.  Ufimtsev  versus  3DT  Resuits:;  Monostatic  Cross  Section  of  a  Wire, 
Circular  Polarization,  2f/\  1.4097 
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Fig.;  1  ’>.  Ufimtsev  versus  SDT  Results:'  Monostatic  Cross  Section  of  a  Wire, 
Circular  Polarization,  2 1/\  S.4113 
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P'ig.  16.  Ufimt3ev  versus  t"  ’T  Results:  Monostatic  Cross  Section  of  a  Wire, 
Linear  Polarization,  2(/\  =  11,0145 
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Fig.  17.  Ufimtsev  versus  SOT  Results 
Monostatie  Cross  Section  of  a 
Wire,  Linear  Polarisation, 


IV.  DISCUSSION 


In  this  paper,  the  theory  of  Van  Vleck,  et  al. ,  has  been  reexamined 
in  some  detail  to  determine  if  the  general  results,  without  the  approxima¬ 
tions,  are  applicable  for  all  angles  of  incidence.  It  has  been  found  that  up 
to  pi  =  157,  the  general  theory  agrees  very  well  with  the  RCS  results  calcu¬ 
lated  by  the  BRACT  computer  program.  At  the  largest  pi  value  considered 
(Pi  =  157),  the  difference  is  less  than  2  dB  at  the  RCS  maxima,  though  it 
is  considerably  larger  at  the  nulls.  The  disagreement  is  smaller  at  the 
maxima  for  thinner  wires.  Calculations  have  been  performed  to  verify 
these  findings.  When  results  obtained  using  the  approximate  and  the  general 
Van  Vleck  theories  are  compared,  it  is  found  that  differences  occur  that 
cannot  be  accounted  for  by  the  neglect  of  angular  dependence  in  the  K,  G', 
G",  H',  and  H"  terms  of  Van  Vleck,  et  al.  These  further  differences  are 
due  to  the  "small"  terms  dropped  from  the  G',  G",  H1 ,  and  expressions 
ant'  not  to  the  use  of  the  asymptotic  expressions  for  the  Cin  and  Si  functions. 
The  correct  expressions  with  these  terms  are 


2G'  = 


+(Pf)(l  -§  F")  -  f  F'E(PI) 
(PI)  +  H*  (pi) 


(36) 


H  (PI)  (l  -  \  Fv)  +  |  F'  tMPD 
4>2  (pi)  +  E  (PI) 


(37) 
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cannot  be  accounted  for  by  the  neglect  of  angular  dependence  in  the  K,  G', 
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2G'  = 


■4>m  ( i  - 1  f")  f's  (pi) 
~~  ‘i'2  (PI)  +  E*  (pi) 


(36) 


H  (pi)  (l  -  f  F--')  +|  F'  «1»(PD 
‘I'2  (PI)  +  E2  (PI) 


(Pi"  2/ 

l-f  F'S( 

Pi  - !) 

>  -  f) 

l  + 

-  I 

) 

(37) 
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